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Abstract. We derive a new shell model of magnetohydrodynamic (MHD) turbulence 
in which the energy transfers are not necessary local. Like the original MHD equations, 
the model conserves the total energy, magnetic helicity, cross-helicity and volume in 
phase space (Liouville's theorem) apart from the effects of external forcing, viscous 
dissipation and magnetic diffusion. In the absence of magnetic field the model 
exhibits a statistically stationary kinetic energy solution with a Kolmogorov spectrum. 
The dynamo action from a seed magnetic field by the turbulent flow and the non 
linear interactions are studied for a wide range of magnetic Prandtl numbers in both 
kinematic and dynamic cases. The non locality of the energy transfers are clearly 
identified. 
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1. Introduction 

Pioneering shell models of hydro dynamic turbulence were developed in the seventies 
[D [21 El IU E], aiming at reproducing the main turbulence features, as intermittency, 
with a low order model of equations. Such shell models were also called wave 
packet representation [6J for the Fourier space is logarithmically divided into shells 
of logarithmic width A such that each wave packet (or shell) k n is defined by k^X 11 ^ 1 < 
k < ko\ n . Those models are local in the sense that each shell interacts with only 
the first neighbours like the DN model (named after Desnyansky and Novikov [1]), 
or the two first neighbours like the GOY model (named after Gledzer, Ohkitani and 
Yamada [31 [7]). The latter has been intensively studied (for a review, see [H El [101 [11] 
and references therein) and subjected to improvements leading to the so-called Sabra 
model [121 [131 E] or extensions using the wavelet decomposition [15] . The GOY and 
subsequently Sabra models have been used in different contexts like convection [16J, 
rotation [T7] or intermittency [TBI [19]. It has been shown [20J that the DN model is 
a spectral reduction of the GOY model, showing in some sense the consistency of one 
model against the other. 

To our knowledge only one non-local shell model of turbulence has been developed 
so far, by Zimin and Hussain [21], projecting the Navier-Stokes equations onto a wavelet 
basis, and reducing the number of variables from statistical assumptions. In such non- 
local model each shell may interact with any other shell. Since then, the original model 
has been improved by Zimin in order to include left- or right-handed polarity of the 
solenoidal basis functions as in the complex helical wave decomposition, and used by 
Melander and Fabijonas [221 [231 123] . 

The extension of shell models to MHD turbulence has been done including 
either first neighbours interactions [251 1261 [271 EE] or two first neighbours interactions 
[291 [301 13T1 [321 133] eventually including Hall effect [3H ES]- However there is a number 
of situations in MHD turbulence in which assuming the locality of energy transfers may 
become somewhat spurious even if the turbulence is considered as isotropic [33]. This 
is true for example when a large scale external magnetic field is imposed leading to 
Alfven waves [36l [371 138] . This problem has been tackled using a MHD shell model 
implementing non local interactions with the externally imposed magnetic field scale 
[28] . However it has been shown recently using an other method [39J that the other 
non local interactions are also important and may rule out the predicted Iroshnikov- 
Kraichnan k~ 3 ^ 2 spectrum. Non local interactions are also at the heart of the dynamo 
problem i.e. when the magnetic field is produced by the turbulent motion instead of 
being externally applied. For example at large value of the magnetic Prandtl number, 
the magnetic spectrum is observed to peak at scales much smaller than the viscous scale 
[4"0~1 |4"T1 W2\ showing a direct evidence of the importance of non local energy transfers. In 
presence of helicity, two possible mechanisms may generate a large scale magnetic field 
as observed in planets and stars: either a local inverse cascade [43] or a non local direct 
transfer from small to large scales as predicted by the mean- field theory [H]. Which 
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mechanism prevails is still not clear. Recently, the importance of non-local interactions 
has been shown in both hydro dynamic [15] and MHD [I6j H7J HH] turbulence. For recent 
reviews on MHD turbulence and the dynamo problem, see e.g. [IHJ ED] 

In the present paper our aim is to introduce a new non local shell model of 
turbulence which can be used either in its hydrodynamic or MHD form. As we 
are interested by MHD applications, we shall introduce the MHD model only, the 
hydrodynamic one being easily deduced from the former, setting the magnetic field 
to zero. Our model can be understood as a non local version of the Sabra model. This 
involves similar rules for complex conjugations and imposes the value of shell spacing 
equal to the golden number. Our first attempt to derive a non local shell model of MHD 
turbulence was based on the Zimin model [21] . However we realized that this model was 
unable to describe the non local interaction between a large scale velocity and two small 
neighboring scales of the magnetic field. We believe that the one which is described here 
is more relevant to actual isotropic MHD turbulence. 



2. General concept 

2.1. The model 

The model is defined by the following set of equations 

U n = ik n [Qn{U, U, a) + Q n (B, B, -a)] - vk 2 n U n + F n , (1) 
B n = ik n [Q n (U, B, b) + Q n (B, U, -&)] - V klB n , (2) 

where 

N 

Qn(X, Y, c) = T rr J\c m X n j rm Y n j rm -\-i + c m X n _ m F n _|-i + c m A n _ m „iY n _i] (3) 

m=l 

represents the non linear transfer rates. The parameters v and i] are respectively the 
kinematic viscosity and the magnetic diffusivity, F n is the forcing of turbulence, and 
k n = X n with A 2 = A + 1 [12J. For iV = 1 in (j3J), we recognize the local Sabra model. 
The additional non local interactions for N > 2 correspond to all other possible triad 
interactions except the ones involving two identical scales. Expressions for the kinetic 
energy E\j and helicity Hjj, magnetic energy Eb and helicity Hb, and cross helicity He 
are given by 

Eu = Y,Mn)> E u {n) = \\U n \\ H v = Y J H u {n) ) H v {n) = \{-l) n k n \U n \ 2 (4) 

n * n * 

E B = Y. E B(n) 1 E B (n) = l -\B n \\ H B = ^2 H B (n), H B (n) = h-l) n k- l \B n \ 2 (5) 

n n * 

H c = Y,H c (n), H c (n) = hu n B: + B n U:)(Q) 

n ^ 

In the inviscid and non- resistive limit (v — r\ — 0), the total energy E = Ejj + Eb, 
magnetic helicity and cross helicity must be conserved (E = H B = H c = 0). This leads 
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to the following expression for the coefficients a l m and b l m : 

1 7 I J. 2 _ ~^m+l ~ ( — I)" 3 _ 



6 1 = (-l) m+1 6 2 = 1 6 3 = -I. (7) 



In the case of pure hydro dynamic turbulence (without magnetic field), the coefficients 
a l m are derived from the kinetic energy and helicity conservations (Ejj = Hjj = 0), 
leading again to the same expression as ([7]). The coefficients T m are free parameters 
depending on m only, that we choose of the form T m = fc^ l _ 1 /A(A + 1). The coefficient 
1/A(A + I) is chosen such that the terms for m = 1 in Q n correspond to the local Sabra 
model. The local Sabra model corresponds to a = — oo. 



2.2. All possible interactions 

In our shell model we see from ([3]) that only a discrete number of triads are allowed. 
For example, Q n does not contain any term involving interactions between the modes 
n + m and n + m + 2. The reason why there is only a discrete number of allowed triads 
comes from the fact that the shells are logarithmically spaced and that the geometrical 
factor A satisfies 

A 2 >A + 1. (8) 

To identify all the allowed triads in a shell model, let us consider three vectors (ki, k 2 , k 3 ) 
satisfying 

k x + k 2 + k 3 = 0. (9) 
Assuming that hi and k 2 belong to shell n and p, we have 

hX"- 1 < |ki| < k X n , koX 13 - 1 < |k 2 | < k X p . (fO) 
From (Q, this implies 

kolX 11 - 1 - A p - X | < |k 3 | < k (X n + X p ). (11) 

Now using the inequalities (JSj) and ( 1TT1) we can show that k 3 belongs to shell q which 
depends on p in the following way: 

p < n — 3 => n — l<q<n+l p = n + 1 =^> q < n + 2 

p = n — 2 =^> n — 2<q<n + l p = n + 2 n<q<n + 3 

p = n — 1 q < n + 1 p > n + 3 p— l<q<p + l 

p = n q < n + 2. 

This is illustrated in Figure [1] in the plane (p, q) where the grey (resp. white) squares 
indicate allowed (resp. not allowed) triads (n,p,q). The demonstration of f|T2|) is given 
in Appendix Appendix A In Figure [H the terms corresponding to the original Sabra 
local model are denoted by "L" , the additional non local terms by "N" and the terms of 
Zimin's model by "Z" . In every case the possible shells (p, q) are symmetric with respect 
to the diagonal in the (p, q) plane . 
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Figure 1. Probability of interactions between three modes belonging to shells n,p 
and q. The white squares correspond to a null probability. 

2.3. Evaluation of a 

There is one free parameter left, a, which corresponds to the non locality strength. It is 
not an easy task (if possible at all) to figure out what a should be in the general case. 
However we tried to estimate it in the case of homogeneous isotropic turbulence (without 
magnetic field). For that, we take two random vectors ki from shell n, and k2 uniformly 
distributed in whole space, and we calculate the probability that k 3 = —(k\ + k 2 ) and 
k 2 belong respectively to shells p and q. It is the simplest estimate of the probability 
that the three modes ki, k 2 and k 3 interact together. A high (resp. low) value of this 
probability is given in Figure [1] by the dark (resp. light) colour of the grey squares. 
In this representation a white square correspond to a null probability. The probability 
spectra along the p and q directions are found to scale as k~ 7 ^ 2 for the " L" and " N" terms 
and £; _5//2 for the "Z" terms (which is consistent with the Zimin's model [21]). In order 
to have terms T m c 2 m and T m c^ in (J3j) scaling in A;" 7 / 2 , we have to take a = —7/2. We 
note that this derivation of a does not imply the terms T m c x m (and neither the diagonal 
terms of Figured]). The latters are determined a posteriori from the conservation laws 
(not included in the probability diagram of Figure [1]) . 

In section 13.1. ll we shall find that the value a = — 5/2 is the most appropriate to describe 
accurately the large scale slope of the kinetic energy spectrum. However, in the rest of 
the paper we shall vary the value of a in order to investigate the role of a in the non 
local interactions. 

2.4- Energy transfers 

To study the non local interactions, we introduce the transfer rate Txy{q, n) from X- 
energy lying in shell q to V-energy lying in shell n. It can operate only within triads, 
implying an additional energy Z lying in all other shells p different from q and n. 
Denoting T(X q \Z p \Y n ) the transfer rate from X q to Y n and which involves Z p as a 
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mediator, the transfer rate from X q to Y n can then be written as 

T XY (q,n) = T{X q \Z p \Y n ). (13) 

From our model (j3J) we see that for each couple (X q , Y n ), only a discrete number of Z p 
can act as mediators. The possible triads are represented schematically in Figure [2j 
Therefore 7xy(g,n) takes the following form 




Figure 2. The three types of possible triad (X q , Z p ,Y n ) : (a) q < n — 1 and p = nil, 
(b) q > n + 1 and p = q ± 1, (c) q = n ± 1 and p < n — 1. Among these cases those 
with repeated subscripts (q = p, p = n or q = n) are forbidden in our model. 



'T(X q \Z n+1 \Y n ) +T(X g \Z n _ 1 \Y n ), for q < n - 2 

T ( ) = \ T ( X l\ Z n+l\ Y n) +Ep<n-2 T ( X q\ Z p\ Y n), fol q = U - I 

xY{q,n) Wpq Zrt+2 | Fn) + E ^_ ir pq Zp | rn)) forg = n + i • { ) 

T(X q \Z q+1 \Y n ) + T(X q \Z q ^\Y n ), for q > n + 2 

Now, we have to define T(X q \Z p \Y n ) for any X q , Z p and Y n . For that, we can re-write 
the model (CQ) (T5]) in the following generic form 

Y n = J2*KM n (X q ,Z p ) + --- (15) 



where M n (X qi Z p ) is a symmetric bilinear form (given in Appendix Appendix B ) 
representing the quadratic terms. The dots represent the dissipation and forcing terms 
appropriate to Y n . 

Now let us denote by S(X q , Z p \Y n ) the rate of energy within the triads (X q , Z p ,Y n ) 
which is transferred from the couple (X q , Z p ) to Y n . It is naturally defined by 

S(X q , Z p \Y n ) = 3? {ik n Y:M n (X q , Z p )} . (16) 

As M n is symmetric, we have 

S(X q , Z p \Y n ) = S(Z p , X q \Y n ). (17) 

In addition, and with the help of (171), we can show that for any triad (X q , Z p ,Y n ) the 
following relation is satisfied 

S(X q , Z p \Y n ) + S(Z P , Y n \X q ) + S(Y n , X q \Z p ) = 0, (18) 

meaning that the energy is conserved within each triad. 

Then we look for T(X q \Z p \Y n ) as a linear combination of S(X q , Z p \Y n ), S(Y n , X q \Z p ) 
and S(Z P , Y n \X q ). As in real MHD turbulence, T(X q \Z p \Y n ) must satisfy the following 
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conditions 



T{X q \Z p \Y n ) = -T(Y n \Z p \X q ) (19) 
S(X q , Z p \Y n ) = T(X q \Z p \Y n ) + T(Z p \X q \Y n ). (20) 



The first condition ( 1191) means that the transfer rate from X q to Y n and from Y n to X q 
with the same mediator Z p are of equal intensity but opposite signs. We note that this 
implies 

T XY (q,n) = -T Y x(n,q). (21) 

The second condition ( 1201) means that in any triad {X q , Z p ,Y n ) the transfer rate from 
{X q , Z p ) to Y n is equal to the sum of the transfer rates from X q to Y n via Z v and from 
Z v to Y n via X q . Then combining (TT7T) . (fT8l) . (1191 and ( 1201) we end up with the following 
expression for T 

TpqZ p |Y n ) = ~ Z p |F n ) - S(y„, Z P |X,)] (22) 

Furthermore we can show that the following energy balances at scale n are satisfied 

Eu(n) + Vu{n) - F(n) = Et 7 "^, n) + T BU (q, n)} (23) 

g 

E B (n) + V B (n) = E Pub (q, n) + T BB (q, n)} (24) 

g 

with the kinetic and magnetic dissipation rates in shell n defined by 

V v {n) = vk 2 n \U n \\ V B {n) = nk 2 n \B n \ 2 (25) 

and 

F(n)= l -{F n U:+F:U n ). (26) 
2. 5. Energy fluxes 

We define the energy flux Iljfy(n) as the rate of loss of X-energy lying in the shells 
j < n to the F-energy lying in the shell j > n. Therefore, we have 

n-l 

n*y(n) = -EEE^(gbli) (27) 
j=o <1 P 
n-l 

= -EEWi)- ( 28 ) 

3=0 9 

The fluxes Ilu^n) and ^^(n) coincide respectively with YJjZl Q{kjU*Qj(U, U, a)} and 
2?=o kjUjQj(B, B, a)}. There is no such coincidence for Ii UB {n) and IT BB ( 



n 



The following flux balances are satisfied 

n-l 

E (^0') + - ^0')) + IW(n) + n B[ /(n) = (29) 

3=0 
n-l 

E (^fl(j) + V B (j)) + TL UB {n) + Il BB (n) = 0. (30) 

j=0 
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In a statistically stationary case they imply that 

n-1 

Uuuin) + IL BU {n) + IL UB {n) + H BB {n) = - £ (V v (j) + V B (j) + . (31) 

3=0 

In order to investigate the importance of non local versus local interactions, we 
define the local part of the fluxes given by ( 1271) in which only the local energy transfer 
rates Txy(n±2|n±l|n), Txy(^±l|^ = Fl|?T-) and Txy(n±l|n±2|n) are involved. These 
fluxes correspond to those of the MHD version of the original (local) Sabra model, taking 
N = 1 in ([3]). The non local parts of the fluxes are defined as the total fluxes minus 
their local parts. 

2. 6. Hydrodynamic forcing 

The hydrodynamic forcing is generally applied at scale n/ and n/+i with nf = 4. It 
is of the form F n = y4 n exp"^ n W where <f> n (t) G [0, 2ir] is constant during time intervals 
t c , the constant value changing randomly from one time interval to the next one. In 
this way we obtain a statistically constant injection rate equal to e = A 2 n t c . We chose 
t c = 10~ 2 for it is smaller than the turn-over time at the injection scale and larger than 
the viscous characteristic time. 

For some arbitrary initial conditions on U n of small intensity we let the 
hydrodynamic evolve until it reaches some statistically stationary state. Then 
introducing at a given time some small intensity of B n we solve the full problem until a 
statistically stationary MHD state is reached. The time of integration needed to obtain 
good statistics depends on v and rj but is typically of several hundreds of time unit (of 
the largest scale n — 0). We define the magnetic Prandtl number as the ratio Pm = u/rj. 

To study the free decaying turbulence, after having reached a statistically stationary 
state, the forcing is set to zero. 

3. Results 

3.1. Energy spectra 

3.1.1. Free decaying turbulence In this section we study free decaying turbulence 
without and with magnetic field. The results are presented in Figure [3] for v = 10~ 6 , 
different values of a and P m . The kinetic (magnetic) spectra are plotted with grey (red) 
dots at different times. The initial conditions are such that the cross-helicity is close 
to zero at any scale. Trying the same simulation but with an initial condition with a 
cross-helicity equal to 1 leads to inertial ranges poorly defined. The time sample at 
which the spectra are plotted is t = 1, 2, 10, 100, 200. The dots corresponding to t = 200 
are the darkest and the smallest. 

We observe that changing a does not change the slope of the inertial range for both 
kinetic and magnetic spectra. This slope compares well with the Kolmogorov slope 
fc~ 2 / 3 (A; -5 / 3 in spectral space) which is represented by the straight line with negative 
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slope in each plot. In the other hand, changing a changes drastically the slope at large 
scale. The kinetic and magnetic energy slopes are indicated by the two straight lines 
with positive slopes in each plot. 

In table 13. lU these slopes are indicated for both spectra and for the different values 
of a that are considered. From hydrodynamic turbulence theory, a slope in k\ (fc 6 in 



a = 


— oo 


-2.5 


-1.5 


-1 


-0.5 


Kinetic slope 


10 


5 


3 


2 


1 


Magnetic slope 


12 


7 


5 


4 


3 



Table 1. Slopes of large scales kinetic and magnetic spectra for several values of a. 

spectral space) is expected at large scales. The only value of a which gives such slope 
is a = —5/2. The slope for the local Sabra model, corresponding here to a = — oo, 
leads to a slope in k^ (k 9 in spectral space) much larger than the one predicted by the 
theory. This is a first drastic difference between the local and non local models which 
emphasizes the importance of including the non local interactions. 

We note that there is always a difference of 2 between the magnetic and kinetic 
slopes, leading to a magnetic spectrum slope in k 7 n (k s in spectral space) for a = —5/2. 
Though a = —5/2 seems to be the most appropriate for hydrodynamic and MHD 
turbulence, in the rest of the paper we shall investigate other values of a as well, in 
order to investigate the role of a in non local transfers. 

3.1.2. Forced turbulence In forced turbulence we consider shells larger than (k n > 1). 
We could also have considered a cut-off at any other arbitrary negative shells. However, 
in order to reach a stationary state it is important to have a cut off scale above which the 
system is not solved (like some integral scale or the scale of the box in an experiment). 
Considering negative shells without limit would imply energy to fill all shells on a time 
scale of the order k~ x and without reaching a stationary state. 

In Figure H] both kinetic and magnetic energy spectra are plotted for three values 
of P m and both kinematic and saturated dynamo regimes. In each plot, the curves cor- 
respond to different values of a. In the kinematic regime (top plots) the Lorentz forces 
corresponding to the term Q n (B,B, —a) are small and can be neglected in the energy 
balance. Then the magnetic energy grows exponentially at any scale. In Figure H] the 
magnetic spectra are normalized by the maximum value of Es{n). In the saturated 
regime (bottom plots), the Lorentz forces act back onto the flow, leading to statistically 
stationary kinetic and magnetic energies. 

In both kinematic and dynamic regimes and for P m < 1 (left and middle plots of 
Figure HJ), the effect of a is not really significant. The kinetic spectra (black curves) are 
almost not sensitive to non local interactions showing that hydrodynamic interactions 
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are mostly local as predicted by the Kolmogorov cascade. Besides the kinetic spectra 
always scale as k~ 2 ^ 3 (fc _5//3 in the spectral space) in the inertial range. We see some 
non local effects onto the magnetic spectra (colour curves) which spread towards large 
scales or even small scales for P m = 10~ 3 . However, the scale for which the spectrum is 
maximum does not change and is roughly equal to k v ~ k u P 3 [ A as again predicted by 
Kolmogorov arguments (see e.g. [33]). 

In the other hand, the non local effect are much more significant for P m > 1 (plots 
on the right), mainly for the magnetic spectra at scales smaller than the viscous scales 
(k n > k u ) and for a > — 1. In the kinematic regime the maximum of the magnetic 
energy spectrum occurs at scales smaller than k v P^[ 4 when a > — 1. In the dynamic 
regime some magnetic bottle neck appears. To understand why it is so, let us first recall 
that the flow scale which produces magnetic field in the most efficient way is the one 
for which the shear is the largest [33]. In the inertial range the flow shear scales as 
k n u n (x k 2 / 3 , and it is then maximum for k n m k u . Therefore the non local interactions 
relevant for the magnetic spectrum extension towards smaller scales are mostly those 
involving U u . The non local terms involving U v and generating magnetic energy Eb(ti) 
with n ^> v involves also P> n ±\. The corresponding non local term in ([2]) are in the 
form k n {T n _ u b 2 n -„U;B n+l + T n _ 1/ -i&^_ l/ „ 1 ?/*-B Tl _i) which scale as Therefore we 

understand that for a + 1 > the non local effect may be strong at small scales. 



P = ID -3 P = 1 P = ID 5 
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Figure 4. Spectra of kinetic (black curves) and magnetic (colour curves) energy 
for three values of P m (from left to right P m — 1CP 3 ; 1; 10 5 ) and several values of 
a (indicated by the labels). The spectra on top (resp. bottom) correspond to the 
kinematic (resp. dynamic) regime. 

3.2. Energy fluxes 

In this section we set a = —5/2 and v = 10 -8 and consider the dynamo saturated 
regime for three values of P m = 10~ 3 , 1 and 10 4 . The kinetic and magnetic spectra 
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are plotted in the top row of figure Both spectra have inertial ranges of Kolmogorov 
types, scaling in k~ 2 ^ 3 (scaling in k~ 5 ^ 3 in the spectral space). For P m = 10~ 3 we identify 
clearly that the magnetic dissipation scale is much smaller than the viscous dissipation 
scale. However the distinction between them is not so clear for P m = 10 4 . 



p = irr 3 p = i p = in 4 




01234567 01234567 01234567 
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Figure 5. Spectra, total and non local fluxes for three values of P m and a = —5/2. 
The kinetic (magnetic) spectra corresponds to black (red) dots. The fluxes Uuu(n) 
and Wbu{ 71 ) ar e represented by blue and green dots. The red dots correspond to 
Huu(n) + ^Bu(n), the black dots to Uuu(n) + Ii B u(n) + U UB (n) + n BS (n). The 
blue, red and black full lines correspond respectively to 5Zj=o ^V(i)> Sj=o ^bU) an( i 

Yr J= o^u{j)+v B (j)). 

The total and non local part of the fluxes are plotted in middle and bottom rows of 
figure [5l The non local part of Ujju{ n ) is found to be always much smaller than Ujju(n), 
implying that the energy transfers are mainly local. In the other hand, the importance 
of the non local part of n S[ /(n) versus the local one depends on P m . 

In figure [6] the ratio il^ nLocaZ (n)/II^ aZ (n) is plotted for the three values of P m . 
For P m = 10~ 3 this ratio is about 20 %. For P m = 1 and for scales smaller than the 
viscous scale k u rs 10 6 , this ratio increases up to 50%. Finally, for P m = 10 4 there is 
a discontinuity at k u , the ratio being then equal to -100% at smaller scales. These are 
first evidences of non local interactions. 
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Figure 6. Ratio U.^ Local (n) /U^ 1 (n) for P m = 1CT 3 (dots), 1 (cross) and 10 4 
(triangle) . 



3.3. Energy transfers 

The four energy transfers are plotted in Appendix Appendix C in figures ICU IC2I and 
IC3I for respectively P m = 10 -3 , 1 and 10 5 . In each figure the column from left to right 
corresponds to a = — oo, — 5/2, — 3/2, — 1 and —1/2. The row from second to bottom 
corresponds to the transfer Tuu(q,n), T BU (q,n), T UB (q,n) and T BB (q,n). The trans- 
fers are plotted versus log 10 q, for three values of n which are indicated by the dashed 
vertical lines on the spectra plots on top row and by the red, green or blue dots. The 
transfers being time-dependent we plot their time-average with error bars corresponding 
to the standard deviation of the mean. This gives some estimation of the robustness 
of the results. Some quantities are much more noisy than the others and then less re- 
liable. The local transfers seem to be always dominant whatever the values of P m or 
a. However there are also some evidence of non local transfers which are discussed below. 



In figure [7] some typical results are presented for the three values of P m = 10 -3 , 1 
and 10 5 (from left to right column), each row from top to bottom corresponding to the 
transfers Tuu(q, n), T B jj(q, n), Tu B (q, n) and T BB (q, n) for one given value of n denoted 
by the black dot. The curves correspond to a = — oo (green), —5/2 (magenta), —3/2 
(blue) and —1 (red). 

• For P m = 10~ 3 , \Tuu\ ^ I^Bfyl? implying that the dominant energy transfer which 
feeds the kinetic energy is a local direct cascade of kinetic energy. 



• For P m = 1 and P m = 10 5 , \7uu\ <^ \T BU \ implying that the kinetic energy is 
mainly obtained from magnetic field. This transfer T B u is found to be mainly local. 



• For P m = 10~ 3 we find that T B y is mainly local and always negative. It means 
that energy is transferred locally from U n to B n (which is consistent with the curve 
T UB just below). In addition, we see that the curve of T BU extends towards larger 
and larger scales when a goes to zero. Though it is small (up to 20%), it is a 
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clear evidence of non local transfer from small scale kinetic to large scale magnetic 
energy. We interpret it as an alpha-effect in the sense of mean field theory. 



• The curves of Tub for P m = 10 5 show also clear evidence of non local transfers. In 
this case, B n is fed by scales of U much larger than n. There is also some local 
transfer back from B n to U n shown by the curves of Tbu- 



• For P m = 10 3 and though much smaller than Tub, there is clear evidence of non 
local direct cascade of magnetic energy as shown by T B b- 



• Finally, for P m = 1 and P m = 10 5 , Tbb shows evidence of some kind of non local 
inverse cascade though much smaller than Tub- 

In figures |DT| ID2I and ID3I of Appendix | Appendix D| we give some qualitative illustration 
of the previous interpretation of the results. 



4. Discussion 



The main originality of the shell model presented in PJPD is that it takes into account 
all possible non local interactions between different shells. In essence it is a non local 
MHD version of the Sabra model. Deriving the model we find that there is one free 
parameter that we call a (and which is different from the alpha-effect of the mean-field 
dynamo theory). In some sense this parameter controls the strength of non locality of 
the model. An estimation of a has been done on the basis of simple probabilistic 
arguments of possible triad interactions, leading to a = —7/2. However, in free 
decaying hydrodynamic turbulence, we find numerically that the value a = —5/2 
gives the right slope of the kinetic spectrum at large scales (fc 4 as predicted by the 
turbulence theory in spectral space). Therefore a = —5/2 seems the most appropriate 
value for hydrodynamic turbulence. In comparison we note that the local Sabra model 
(corresponding here to a = — oo) leads to a slope about k 9 . 

In the other hand, in MHD turbulence it is not obvious that a has to be the same 
as in hydrodynamic turbulence. Actually, it is even not obvious that a should not 
depend from the scale n, which is not possible in our model. Therefore we investigated 
several values of a ranging several values from — oo to zero. Several values of P m have 
been investigated depending if it is lower than, equal to or larger than 1. In order to 
characterize the non local energy transfers from X q to Y n within any possible triad, 
the quantities T X y(q,n) have been derived an plotted versus q for a few values of n. 
Though most of them are local, several energy transfers have been found to be partially 
non local, depending on P m . 
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Figure 7. From left to right column, P m = 10 -3 , 1,10 s . From top to bottom, the 
plots correspond to Tjju(q,n); Tbu(q, n); Tub(q, n); Tbb{<1, n) for n indicated by the 
black dot. Each curve correspond to a = — oo (green), —5/2 (magenta), —3/2 (blue) 
and —1 (red). 



Acknowledgments 

This work was supported by the Econet Program 10257QL. Multiprocessor 
computational resources were supported by special program of the Russian Academy of 
Science. RS thanks BRHE post-doctoral fellowship program (CRDF grant Y2-MP-09- 
02) for financial support. The authors thank A. Alexakis, P. Frick, H. Politano and A. 
Schekochihin for fruitful discussions. 

Appendix A. Possible triads in logarithmic shell models 

For a given shell n, we first take benefit from the fact that the triads (n, p, q) and (n, q, p) 
are identical and then the representation in the plane (p, q) is symmetric with respect 



A non local shell model of turbulent dynamo 



16 



to the diagonal. We then limit our demonstration to half of (TT21 : 

• Assuming p < n — 1, we have A" + X p < X n + A™" 1 . Then from ([8]) and (TTTj) we 
have | ki3 1 < A;oA n+1 , implying q < n + 1. 

• For p = n, we have A" + A p = 2 A™. From ([8]) we have 2 < A 2 . Therefore from ( TTTT) 
we have |k 3 | < /c ^ n+2 an d then q < n + 2. 

• From (ED we have A p+1 - X p > X p ~\ For p > n + 1 we then have A p+1 - A p > A n , 
implying from (TXT]) that |k 3 | < A; A P+1 and then q < p + 1. 

• Assuming that p < n — 3, from (TTTj) we have |k 3 | > /c |A n_1 — A ra_4 |. Then (JH]) 
implies A™" 1 > A n " 2 + A n ~ 3 > A™" 2 + A n " 4 . Then we have |k 3 | > k X n ~ 2 > A™" 1 
and then q > n — 1. 

• For p = n — 2, from ( TTTT) we have |k 3 | > fco|A n_1 — A n ~ 3 |. From ([8]) we have 
A n_1 — A n_3 > A"~ 2 , implying |k 3 | > A™ -2 and then q > n — 2. 

Appendix B. Expressions of the energy transfers 

The energy transfer ( |22l) can be written with the help of ( flBl) in the form 

T(X q \Z p \Y n ) = ~3i {ik n Y*M n {X q: Z p ) - ik q X* q M q (Y n , Z p )} . (B.l) 
The symmetric bilinear form M n (X q , Z p ) is defined as follows 

M n (U q , U p ) = L n {U q , U p , +a) + L n {U p , U q , +a) (B.2) 

M n (B q , B p ) = L n {B g , Bp, -a) + L n (B p , B q , -a) (B.3) 

M n (U q , Bp) = L n (U q , B P1 +b) + L n (B p , U q , -b) (B.4) 

M n (B q , Up) = L n (B q , U p , -b) + L n (U p , B q , +b) (B.5) 

with L n (X p , Y q , c) given by 

{T p _ q c p _ q X q Z p for q < n — 2 and p = n — 1 

T n - q cl_ q X*Z p for q < n - 1 and p = n + 1 .(B. 

T q _ n c q _ Tl X*Z p for g > n + 1 and p = g + 1 

We note that L n (X p , Y q , c) is related to Q n (X, Y, c), defined in ([3]), by 

Q n (X,Y,c) = J2 L n(X P ,Y qj c). (B.7) 

Appendix C. Energy transfers results 
Appendix D. Illustration of the energy transfers 
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Figure Dl. Illustration of the energy transfers for P m <C 1 for local (top) and non 
local (bottom) model. The thickness of the arrow gives some qualitative estimate of 
the strength of the transfer. 
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